JOURNAL OF RESEARCH of the National Bureau of Standards— A. Physics and Chemistry 
Vol. 68A, No. 6, November-December 1964 

Energy Dependence of the D-D Reaction Cross Section 

at Low Energies ^ 

J. G. Brennan^ and J. J. Coyne 

(August 11, 1964) 

The energy variation of th(^ approach cross section for the D-D reaction is calculated 
using several different methods. The simplest method assumes strong absorption inside 
the nuclear surface and uses the WKB approximation. The slope of the Gamow plot is in 
this case a constant. A strong absorption model, without WKB approximation, leads to 
negative corrections to this constant slope. A weak absorption model, which seems more 
applicable, gives both positive and negative corrections, depending on the depth of the 
nuclear potential. Finally, sinc(^ the amount of absorption seems important, the effects of 
varying the absorption are studied using the optical model. The possible influence of the 
P- waves is also investigated. An experimental study of the energy dependence of the slope 
of the Gamow plot should give information on th(^ optical potential which one deuteron sees 
when it approaches another deuteron. 



1. Introduction 

In an earlier treatment [1] ^ of the D-D reaction 
cross section, the strong energy dependence at low 
energies was attributed to the effects of the Coulomb 
and centrifugal barriers on the various angular 
momentiHu components of the incident beam. How- 
ever, the barrier penetration factors were calculated 
using the WKB approximation for the (voulomb 
wave functions and retaining only the main energy 
dependence. If, in addition to these approxima- 
tions, we also go to the low energy limit, then the 
energy dependence of the reaction cross section is 
given by the simple expression 



an={AIE) exp i-BE-^) 



(1) 



where A and B are hidependent of the energy. A 
plot of the logarithm of (TrE as a function of E~^ 
would give a straight line under these approxima- 
tions. Such a plot is often called a Gamow plot 
and is frequently used in the analysis of the low 
energy data. 

Looking at the experimental results from three 
groups [2-4] in the energy region of about 8 to 100 
kev (center of mass) , Brennan [5] has indicated that 
the data might suggest an increase in B of eq (1) 
with energy. Brennan calculated the D-D reaction 
cross section using a low energy approximation [6] 
to the small absorption method of Ostrofsky, Breit, 
and Johnson [7]. In this present paper, the energy 
dependence of the reaction cross section has been 
studied under various assumptions suggested by this 
previous work. 
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wSince retaining ojily the main energy dependence 
in the WKB a])])roximation is equivalent to assimi- 
ing the behavior of the wave function due to scatter- 
ing from inside the nuclear surface does not affect 
the energy dependence of the cross section, the 
reaction cross section was first calculated using the 
continuum model [(S], which assumes complete absorp- 
tion inside the nucleus and, tlierefore, does not allow 
coherent scattering from inside the nuclear surface. 
The results of tliis calculation are compared with the 
results of reference 1. Th(Mi])])roach used by Bren- 
nan, based on the small absorption model, involves 
an approximation valid only at low energies. To in- 
vestigate the validity of this approximation, a sec- 
ond calculation was done, in which this approxima- 
tion was not made. 

Since these two calculations involved models 
assuming strong and weak absorption, the effects 
of absorption were studied in a third calculation, 
using the optical model [9]. An attempt was also 
made to study the influence of P- waves on the cross 
section although only S-waves were included in the 
three main calculations mentioned above. 

2. General Considerations 

In this paper, we are considering the D-D reaction 
cross section at energies below 1 MeV in the center- 
of-mass coordinate system. In this low energy 
region, it is valid to treat the initial kinetic energy 
as small compared with the energy of the nucleons 
once they are inside the nucleus. Since the energy 
dependence of the reaction cross section does not 
show any of the sharp peaks which are found when 
there is compound nuclear resonance, we can also 
assume that there are no resonances to compound 
nuclear states below 1 MeV. 

Under these two assumptions, it is a good approxi- 
mation to assume that the reaction cross section 
factors into two parts. The first part, the approach 
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cross section, contains all of the energy dependence 
and measures the probability of the bombarding 
particle getting inside the target mucleus. The 
second part, the specific nuclear factor, is assumed to 
be independent of the energy and is related to the 
probability that, once the bombarding particle gets 
inside the target nucleus, the system will make a 
transition from the initial to the final state. 

The reaction cross section can also be subdivided 
into cross sections for reactions initiated by particles 
with different I values, i.e., 



<^/e = S(4/9)Pz(Tz 



(2) 



where Pi is the specific nuclear factor for initial 
orbital angular momentum / and ai is the corre- 
sponding approach cross section. The factor 4/9 
arises from considerations of spin and symmetry [1] . 
Following previous treatments [1], we have ignored 
the contributions from the quintet spin states. In 
the present treatment, no attempt has been made to 
calculate the specific nuclear factors, Pi. The}^ are 
assumed independent of the energy, and, therefore, 
do not affect the energy dependence of aR. 

3. WKB Approximation and the Continuum 
Model 

The continuum model [8] assumes complete 
absorption inside the nuclear surface and thus 
eliminates the effect of coherent scattering from 
inside the nucleus. The problem of calculating the 
approach cross section becomes one of calculating a 
potential barrier transmission coefficient. The ap- 
proach cross section is related to the transmission 
coefficient as follows : 



ai=ir^\2l+\)Ti 



(3) 



where 7rX^(2/+l) is the geometric limit and Ti is the 
transmission coefficient. The transmission coeffi- 
cients are given in terms of A ^ and /S/, the penetration 
and shift factor, by Blatt and Weisskopf [10]. 
A^ and Si are, in turn, expressed in terms of the 
regular and irregular Coulomb wave functions and 
their derivatives [11]. 

If the WKB approximation is used for the Coulomb 
wave functions, and only the main energy dependence 
is retained, then we have [1] 

(Xi=Qx^[—2Ci] 
Ci=r\ki{r)\dr', 

JR 

ki{r)={2ixl\AE-e'lr-h'{l + hyi2y.r']Y^ (4) 

and Vi is the value of r at which ki{r) = {). It should 
be noted that in excluding all but the most important 
energy dependence, we have also eliminated all 
dependence on TJq^ the depth of the nuclear potential] 



The approach cross section, calculated using the 
WKB method, depends only on R, the nuclear 
radius, which we have chosen to be 7X10"^^ cm 
throughout this paper. This value of R is the one 
commonly used in the D-D reaction calculations. 
Unless R is changed drastically, the only effect of 
varying R would be to change the height of the 
potential barrier. 

Two calculations have been performed for the 
S-wave approach cross sections, one using the WKB 
approximation given by eq (4), the other using exact 
Coulomb wave functions [12] in the continuum model. 
In the second calculation, which we shall refer to 
simply as the continuum model, the dependence on 
Uq was also retained. Figure 1 gives the slope of 
the Gamow plot for the WKB calculation and for 
the continuum model calculations with various well 
depths. In this figure and in all subsequent graphs 
of the slope of the Gamow plot, Ec corresponds to 
the height of the Coulomb barrier, 205.69 keV; 
the scale marks a and h corresponds to 100 keV and 
8 keV, respectively. The depths of the nuclear 
potentials are given on the curves in MeV. The 
results of figure 1 were obtained using numerical 
differentiation and using double interpolation of 
tabular data [12], good only to four significant 
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Figure 1, B(E) versus E~i, WKB approximation and the 
continuum model. 

In this and all subsequent figures, the numbers on the curves are the depths of 
the nuclear potential in MeV, Ec corresponds to the height of the Coulomb 
barrier, a corresponds to 100 keV and 6 to 8 keV. 



676 



figures. Consequently the accuracy of the curves 
shown is only about 10 percent. Except for Uq=0, 
which is certainly an unrealistic value for the depth 
of the nuclear potential, all of the continuum model 
curves give negative corrections to the WKB results. 
In order to test the importance of higher angular 
momentum, the WKB calculations were performed 
also for /= 1, 2, and 3. If we assume all /-waves have 
the same specific nuclear factors, then, for compari- 
son, we can simply add the cross sections. Figure 2b 
shows the effect of adding in higher and higher / 
values, while figure 2a compares the slope of the 
Gamow plot for do with that for (0-0+ (Ti). Higher 
Z- waves were not included in figure 2a because the 
results were not significantly different in the region 
of interest. While the addition of the P-waves can 
cause a positive correction to the slope of the Gamow 
plot, this correction is very small compared to the 
correction indicated by experiment. Moreover, in 
a later and a more accurate calculation involving 
S- and P-waves, we shall see that the relative iiiipor- 
tance of higlicr / values will depend very mucli on 
the depth of the nuclear potential. 




Figure 2a. B(E) versus (E)-2 for S waves and for (S + P) 
waves in the WKB approximation. 



Figure 2b. So-^ versus E in the WKB approximation y showing 
the effects of higher I waves. 



4. Small Absorption Approximation 

The treatment outlined in tlie previous section 
implied complete absorption inside the nucleus, 
since it ignored the effect of colicrent scattering 
from inside the nucleus. The metliod described in 
this section assumes that the absorption is small, in 
fact so small that its effect on the probability density 
inside the nucleus can be ignored. This j net hod was 
first introduced by Ostrofsky, Breit, and Johnson [7] 
and shall be referred to as the OBJ method. We 
will assume the absorption probability is the same 
throughout the nucleus, so that the sliape factoi' for 
absorption will be one. This, we shall see latei*, is 
equivalent to choosing the imaginary part of the 
optical potential to be a square well. The approacii 
cross section is the integral of the probability density 
over the nuclear volume divided by the flux of the 
incident beam. Since, in this method, the approach 
cross section is really a measure of the probability of 
being inside the nucleus, the specific nuclear factor 
is not a ^'branching ratio" but is related to the 
absolute probabihty of absorption. 

The wave function inside the nucleus vvdll be a 
solution for tlie nuclear potential which has been 
assumed. Outside tlie nucleus, the wave function 
is a linear combination of regular and irregular 
(\)ulomb wave functions, chosen so that the wave 
function and its derivative are continuous at the 
nuclear surface. With the wave function determined 
in this way, the approach cross section is given by 
eq (6) of Ostrofsky, Breit, and Johnson [7J. Using 
the OBJ method, we have treated only the S-waves, 
using a square well potential and also a clipped 
harmonic oscillator potential. The Coulomb wave 
functions foi* })ot}i calculations w^ere the same as 
those used for the 'Continuum modeP' calculation. 
The square well calculation was carried out on a 
Burroughs El 01 at the Georgetown College Ob- 
servatory. The calculation for the harmonic oscil- 
lator was done on an IBM 704, located at the U.S. 
Naval Ordnance Laboratory. 

Figure 3 gives Gamow plots for various depths of 
the nuclear potential. Besides the square well and 
the oscillator results, we have also included on each 
graph a plot of the WKB results normalized so that 
they correspond to the square well results at low 
energies and zero depth of the nuclear potential. 
The WKB results are included so that the variation 
of the curves with the depth of the nuclear potential 
will be more obvious. 

While the results for a given depth may be quite 
different for the two potentials, the general behavior 
as a function of Uq is the same. At a given depth 
the square well will give a curve of a certain shape, 
which can be duplicated by a harmonic oscillator 
well of greater depth. As we might suspect, the 
energy dependence of the low energy reaction cross 
section will not give us information of the shape of 
the nuclear potential, but will only indicate what 
depth must be chosen with a given shape, in order to 
get the same results. 
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Figure 3. Gamow plots in the OBJ method for a square well 
and for a clipped harmonic oscillator nuclear potential, for 
various values for the depth, Uq. 



Figure 4 is a graph of the slope of the Gamow plot 
as a function of E~^ for a square well of various 
depths. The WKB results are again included for 
comparison. This graph should be compared with 
figure 1 in order to compare the results of the 
^^continuum'^ model and the OBJ method. The 
comments about the accuracy of the curves in figure 
1 apply equally well to figure 4. However, it is still 
clear that the OBJ method allows both positive and 
negative corrections to the WKB result, while the 
continuum model, for reasonable values of Uq, 
allows only a negative correction. 

5. Small Absorption Approximation at Very 
Low Energies 

Johnson and Jones [6] have shown that in the low 
energy limit the approach cross sections of the OBJ 
method are given by 



<Ti^(guE)(e'-v-l)-\l+U/v') 



(5) 



where rj = e^lhv(xE~^ and qi and ti are energy in- 
Mependent quantities, which are functions of the 
depth of the square well potential. 




OBJ METHOD 
WKB METHOD 



Figure 4. B(E) versus E-| for a square well in the OBJ 
method. 



The expression above was that applied to the 
D-D reaction cross section by Brennan [5]. He 
assumed that he was working at low enough energies 
that only the S-waves were important. If we con- 
sider only S-waves, the slope of the Gamow plot is 
given by 

d(lnaE)ldE-y^= - [(2m) -^2xh/a] (1 +tofTv') (6) 

which corresponds to eq (5) of Brennan [5]. The 
right side of eq (6) is just the slope of the Gamow 
plot, B. The constant corresponds to the WKB 
result, while the second term includes a correction 
proportional to E^^^. In figure 5 we plot ^o a.s a 
function of Ul (the other two curves will be referred 
to later). We see that the correction term can vary 
quite a bit depending on the value of C/q- Analyzing 
the available data [2-4], Brennan [5] decided on a 
value of Uq in the vicinity of 9 MeV. 

In order to show comparative results for the OBJ 
method, and the low energy limit, figure 6 gives the 
cross sections for both at different values of the 
depth. If we compare the values of B calcualted, 
we find that, in general, the low energy limit (6) is 
valid. However, when ^o is large, the results can 
differ considerably even at fairly low energies. 
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Figure 5. to and t of equations (6) and (8) versus U02. 

In order to see how important the P-waves could 
be, the low energy Innit (5) was also applied to the 
P-wave approach cross section. The results, how- 
ever, are somewhat dependent on the ratio of the 
specific nuclear factors for the S- and P-waves. 
Including the P-waves, the total reaction cross sec- 
tion becomes 



(rR= 



4Po2o 



9E(e'-^-l) 



['+l+ftlO+?)]- <'> 



In terms of a Gamow plot this gives 



d{\ncE)/dE-''= - [ (2^) -^2,rh/a] (1 + t/wv^) (8) 



where 



<o+(Pigi/iPogo)<i 
" l + (Px2./Po2o) 



(9) 



The quantity of interest is t and as we see, it is a 
function of jPi/Pq. Figure 5 gives t as a function of 
Ul for Pi/Po=l.O and 0.5. Comparing t to to, we 
see that the P-wave contribution can be very im- 
portant or very unimportant, depending on the depth 
of the nuclear potential. In the vicinity of 9 MeV 
the P-wave is not very significant. 
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Figure 6. o-q versus E for the OBJ method and for the low 
energy limit, for various values of Uq. 



6. Optical Model 

The methods for calculating the approach cross 
sections, which we have treated thus far, are at two 
extremes. Either we assume that the absorption 
inside the nucleus is complete or we assume tliat the 
absorption is so small that it is negligible compared 
to the scattering. In order to study the effect of 
varying the absorption, the S-wave approach cross 
sections were calculated, using the optical model [9]. 
The optical model leads to an absorption which is 
proportional to the probability density inside the 
nuclear surface. Since, as stated previously, Pp, 
the specific nuclear factor of the OBJ method, is 
proportional to the absolute probability of absorp- 
tion, and Vcy the imaginary part of the optical 
potential, is proportional to the absorption in the 
optical model, the two methods sliould give the 
same results for small absorption. Comparing the 
reaction cross section for the optical model with 
that for the OBJ method, we find that, in the limit 
of small absorption, 



Po=^2VJw. 



(10) 
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Figure 7. B(E) versus (E)-i, the optical model, f, the 
ratio of the imaginary to the real part of the potential= 10. 



Calculations in the optical model were carried out 
on the IBM 704 at the National Bureau of Standards. 
The potential used was a square well of the form 
Uo{l+i^). However, the Coulomb wave functions 
used were calculated directly by the computer and, 
therefore, are probably better than the Coulomb 
wave functions used in the previous calculations. 
S-wave approach cross sections were calculated for 
various values of Uq and f. Figures 7 through 11 
give B, the slope of the Gamow plot as a function of 
E~K Comparing figure 11 with figure 4, we see that 
the results for f= 0.0001 and the OBJ method com- 
pare as favorably as we can expect. Calculations 
were also done for f=0.0002, but the cross section 
results differed from f= 0.0001 results by only a 
factor of 2, while the B results coincided in all cases. 
This indicates that when the absorption is low 
enough to make r= 0.0001, the OBJ method is 
surely a good approximation, i.e., the absorption 
does not affect the approach cross section. 

There is no continuum model equivalent in the 
optical model, for as the imaginary part of the poten- 
tial increases, the absorption inside the nucleus 
increases, but the reflection at the nuclear surface 



also increases. As the imaginary part of the optical 
potential approaches very large values the reaction 
cross section goes to zero. However, a comparison 
of figures 1 and 7 shows that the shape of the B 
versus E~^ curves is roughly the same for the continum 
model and the optical model with large values for f . 



7. Summary and Conclusions 

The experimental results for the D-D reaction 
cross section at low energies are usually given in 
terms of a straight line on a Gamow plot. Using 
the WKB method, we have found that the calculated 
value for B(E), the slope of the Gamow plot, is 
indeed very close to constant in the energy region 
from 8 to 100 keV. The ' 'continuum'' model, which 
is the more accurate strong absorption approxima- 
tion, gives a negative energy dependent correction to 
the constant value predicted for B(E) by the WKB 
method. Brennan [5] has pointed out that the 
experimental results seem to indicate a positive 
correction to B{E). In addition, theoretical con- 
siderations seem to indicate that a weak, rather than 
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Figure 8. B(E) versus E~2, optical model, ^=1. 
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Figure 9. B(E) versus E— 2-, optical rnodel, ^=10-^. 



a strono-, absorption model is more applicable to the 
D-D reaction. Using- the OBJ method, wliicli 
ignores absorption in calculating the approach cross 
section, we have found that botli positive and nega- 
tive corrections to B{E) are possible. (\)mparing 
the results for a square well and a clipped harmonic 
oscillator imclear potential, we have found that the 
approach cross section is not very sensitive to the 
shape of the nuclear potential. 

To study the effects of absorption we have also 
calculated approach cross sections using the optical 
model. When the imaginary part of the optical 
potential becomes small, the optical model results 
approach the results of the OBJ method, as expected. 
There is no limit of the optical model which cor- 
responds to the ^'continuum" model, but as the 
imaginary part of the optical potential becomes 
large, the shape of the results approach those of the 
^'continuum" model. 

The calculations discussed so far assumed that only 
the S-waves were important. Using a low energy 
approxhnation to the OBJ method, a study was made 
of the possible influence of P-waves. It was found 




E,( kev)^ 



Figure 10. B(E) versus E-§, optical model, f=10-2. 



that tlie hnportance of the P-waves was very much 
a function of the depth chosen for the potential. 
More accurate experimental results for the total 
reaction cross section will, we hope, lead to a better 
knowledge of the energy dependence of B(E), which 
in turn would tell us more about the optical potential 
proper to the collision between two deuterons. 
Also, more accurate studies of the energ}^ dependence 
of the angular distribution of the reaction products 
might be used to investigate the influence and im- 
portance of the P-waves. 



The authors express their gratitude to the U.S. 
Naval Ordnance Ijaboratory, White Oak, Silver 
Spring, Md., to the Georgetown College Observatory, 
Washington, D.C\, and to the National Bureau of 
Standards, Washington, D.C., for the generous use 
of their computer facilities. 



681 
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Figure 11. B(E) versus E~i, optical model, ^=10~^. 



8. References 

[1] F. M. Beiduk, J. R. Pruett, and E. J. Konopinski, Phys. 

Rev. 77, 622 (1950). 
[2] Arnold, Phillips, Sawyer, Stovall, and Tuck, Phys. Rev. 

93, 483 (1954). 
[3] Eliot, Roaf, and Shaw, Proc. Roy. Soc. (London) A216, 

57 (1953). 
[4] Davenport, Jeffries, Owen, Price, and Roaf, Proc. Roy. 

Soc. (London) A2U, 66 (1953). 
[5] J. G. Brennan, Phys. Rev. Ill, 1592 (1958). 
[6] J. L. Johnson and H. M. Jones, Phys. Rev. 93, 1286 

(1954). 
[7] M. Ostrofsky, G. Breit, and D. P. Johnson, Phys. Rev. 

49,22 (1936). 
[8] H. Feshbach, D. C. Peaslee, and V. F. Weisskopf, Phys. 

Rev. 71, 145 (1947). 
[9] H. Feshbach, C. E. Porter, and V. F. Weisskopf, Phys. 

Rev. 96,448 (1954). 
[10] J. M. Blatt and V. F. Weisskopf, Theoretical Nuclear 

Physics, p. 360 (John Wiley and Sons, Inc., New 

York, N.Y., 1957). 
[11] Ibid., p. 332. 
[12] I. Bloch, M. M. Hull, A. A. Broyles, W. G. Bouricius, 

B. E. Freeman, and G. Breit, Rev. Mod. Phys. 33, 

147 (1951). 



(Paper 68A6-317) 



682 



